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ABSTRACT 

Wind-interaction models for gamma-ray burst afterglows predict that the 
optical emission from the reverse shock drops below that from the forward shock 
within 100s of seconds of the burst. The typical frequency v m of the synchrotron 
emission from the forward shock passes through the optical band typically on a 
timescale of minutes to hours. Before the passage of v m , the optical flux evolves 
as t -1 / 4 and after the passage, the decay steepens to £-( 3 p- 2 )/ 4 ; where p is the 
exponent for the assumed power-law energy distribution of nonthermal electrons 
and is typically ~ 2. The steepening in the slope of temporal decay should be 
readily identifiable in the early afterglow light curves. We propose that such 
a steepening was observed in the R-band light curve of GRB 021004 around 
day 0.1. Available data at several radio frequencies are consistent with this 
interpretation, as are the X-ray observations around day 1. The early evolution 
of GRB 021004 contrasts with that of GRB 990123, which can be described by 
emission from interaction with a constant density medium. 

Subject headings: gamma-rays: bursts — stars: mass loss — stars: supernovae: 
general 
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1. INTRODUCTION 

The initial model for the afterglows of gamma-ray bursts (GRBs) involved synchrotron 
emission from nonthermal electrons accelerated to a power-law spectrum in a relativistic 
spherical blast wave expanding into a constant-density, presumably interstellar, medium 
(ISM; Meszaros & Rees 1997). This simplest model has difficulty, however, explaining 
quantitatively the dozen or so sources whose afterglows are observed well enough to allow 
for detailed modeling (e.g., Panaitescu & Kumar 2002). The most commonly discussed 
complication is a collimated energy injection (Rhoads 1997; Sari, Piran & Halpern 1999). 
This jet model provides a reasonable fit to the majority of the well observed afterglows, 
assuming a const ant- density ambient medium (Panaitescu & Kumar 2002). For some 
sources, a blast wave expanding into an ambient medium of r~ 2 density distribution, as 
expected of a stellar wind, can fit the data equally well or even better (Chevalier & Li 
2000, CL00 hereafter; Li & Chevalier 2001; Panaitescu & Kumar 2002). Notable examples 
include GRB 970508 (see Frail, Waxman & Kulkarni 2000 for a different view) and GRB 
011121 (Price et al. 2002). There is evidence for two types of GRB ambient environments, 
with implications for their progenitors. They are not immediately distinguishable because, 
at an age of a few days, the preshock wind density is comparable to an interstellar density. 
At earlier times, the density contrast is higher and the jet effects are less important. The 
early afterglow observations are expected to provide a better handle on the nature of the 
ambient medium. 

In this Letter, we summarize the characteristics of the early afterglows expected in the 
wind model (§ 2), and argue that the afterglow observations of GRB 021004 are consistent 
with the source interacting with a Wolf-Rayet type wind (§ 3). The strongest evidence 
for wind interaction comes from the initial slow decay of the R-band light curve and its 
prominent steepening around day 0.1. 
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2. ANALYTIC LIGHT CURVES OF EARLY AFTERGLOWS 

Analytic light curves for the standard ISM model are given in Sari, Piran & Narayan 
(1998), assuming a power-law electron energy spectrum with index p and constant 
fractions of blast wave energy, e e and e#, going into nonthermal electrons and magnetic 
field respectively. At any given time, the synchrotron spectrum is determined by a 
set of characteristic frequencies: the typical frequency v m , cooling frequency v c , and 
self- absorption frequency v a . The light curve at any given frequency v is determined by the 
characteristic times t m , t c and t a , when v m , v c and u a cross v, and the critical time t when 
v m and v c become equal. The light curve of the ISM model was extended to the wind case 
by CLOO (see also Panaitescu & Kumar 2000 and Granot & Sari 2002). In the wind model, 
the cooling and self-absorption frequencies are expected to be lower at early times than 
those in the ISM model, because of a higher ambient density at small radii. The cooling 
frequency has a strong effect on the early emission in the optical and the self-absorption 
frequency in radio. 

At the earliest times, the optical flux may be dominated by the emission from the 
reverse shock front, but the forward shock front is not much fainter for wind interaction 
(CLOO, eq. [58]), assuming that the two shocks have similar efficiencies for the production 
of synchrotron radiation. Once the reverse shock front has passed through the initial 
shell, which is expected to occur on a timescale similar to that of the gamma-ray burst, 
the reverse shock emission drops sharply because it is in the fast cooling regime. Unless 
the reverse shock is "refreshed" by a continued flow, the decline is determined by off-axis 
emission that arrives at a later time (Kumar & Panaitescu 2000). Optical frequencies are 
expected to be between v c and u m , so the flux is ~ z/^ 1 / 2 and decays as t _5//2 (Kumar & 
Panaitescu 2000). 

To describe the forward shock emission, we use the characteristic times mentioned 



- 5 - 



above, which were estimated by CLOO assuming p = 2.5. We rescale the estimates to the 
R-band with v R = (z//4.5 x 10 14 Hz), and obtain 

t m = 0.04(1 + zf'hf^e^El^- 2 '" days, (1) 

which is about an hour for standard parameters. The parameter z is the cosmological 
redshift, E 52 is the blast wave energy in units of 10 52 ergs, and e n = e/10™. Note that t m 
does not depend on the wind density A directly, where p = Ar~ 2 , although the wind must 
be sufficiently dense that the transition from fast cooling to slow cooling occurs after t m . It 
depends most sensitively on e e . The condition on e e for the R-band break to occur between 
one minute and one hour is 0.005 < e e (l + zf^e^Elj^ < 0.1, which covers a reasonable 
range, and has only a weak dependence on other parameters. The transition to slow cooling 
occurs around the time to — 1(1 + z)e e -\eB,-\A* day, which is typically later than t m in the 
R-band. Here, A, = A/5 x 10 11 g cm -1 . The cooling frequency crosses the frequency u R at 
a time t c = 5 x 10 3 (1 + 2) 3 e^ _ 1 i? 52 1 A 4 z/ 2 ? days, which is typically later than both t m and 
to- The self- absorption frequency u a is typically well below the optical, and can be ignored. 

Before t , the synchrotron electrons are in the fast cooling regime, and the flux peaks 
at the cooling frequency u c , so F Vc = F u>max with 

F^ max = 2.1(1 + z^e^E^A^t'^ mJy, (2) 

where dn is the luminosity distance in units of 10 Gpc. The flux at the typical frequency 
v m is lower, and is given by 

F Vm = 2.7(1 + z^X^lEU'dlft^ mJy, (3) 

which is independent of the wind density. After the transition to slow cooling 
at t , the flux peaks at v m rather than u c , and now F Vm = F vmax (eq. [2]). 
There is a general scaling for the late-time R-band flux after the break (t > t m ) 

F UR (t) OC er l e g- 2 )/ 4 ^+ 2 )/ 4 t -(3 P -2)/4 K i 3(p-l)/4 e -l/4 E 3/4 r( 3 p _ 2)/ 4 j where ^ ^ ^ break 

time at the R-band. Sources with later breaks tend to be brighter. 
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3. MODELING GRB 021004 

The gamma-ray burst GRB 021004 was detected by the HETE II satellite (Shirasaki et 
al. 2002) and had an optical afterglow detected within minutes of the 7-ray burst (Fox et 
al. 2003b). The early light curves of the only other afterglows detected at such early times, 
GRBs 990123 (Akerlof et al. 1999) and 021211 (Fox et al. 2003a), have similar shapes, 
both showing a rapid initial decline followed by a slower decay, although GRB 021211 
is fainter by about 3 magnitudes at similar epochs. In both sources, the initial rapidly 
decaying emission is interpreted as coming from the reverse shock of GRB ejecta running 
into a constant density medium (Sari & Piran 1999; Li et al. 2003; Fox et al. 2003a; Wei 
2003). The early optical afterglow of GRB 021004 shows a different behavior: it has a very 
slow initial decay of t~ 0A±0A , followed by a steepening around day 0.1 into approximately 
t L (Fox et al. 2003b). Kobayashi & Zhang (2003) interpreted the early afterglow data 
in terms of emission from a combination of reverse and forward shocks expanding into a 
constant-density medium. Fox et al. (2003b) questioned this interpretation, and suggested 
instead a continued energy injection into the blast wave after the 7-ray burst to explain 
the initial slow decay. We propose that the slow decay is a natural consequence of the fast 
cooling [y c < u m ) expected in a wind model at early times and that the steepening is caused 
by the typical frequency u m passing through the optical band from above while still in 
the fast cooling regime. We show that this interpretation, besides fitting the R-band light 
curve, is in a reasonable agreement with the radio data available at several frequencies. 

The free parameters that appear in the wind model can be estimated analytically 
using the formulae given in the previous section. Our identification of the break 
time t m ^ = 0.1 day in the R-band light curve yields, using eq. (1), the relation 
e ^-i4 /3 -i^52 3 = 1-7, for a redshift of z = 2.32 (Chornock & Filippenko 2002). For a 
cosmological model with H = 71 km s _1 Mpc -1 , tt M = 0.27, and f2 A = 0.73, this redshift 
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corresponds to d L1 = 1.89. At the break, the R-band flux is approximately F Vm = 0.83 mJy, 
which yields a second relation ^-x^-x^a =1.9 from eq. (3). The wind density A* does 
not enter into either of the two relations, which enables us to express e e and E 52 in terms of 
e B : e e = O.lle^ and E b2 = 4.0e^ 3 _ 1 . 

To constrain the wind density A*, we note that the typical frequency v m decreases with 
time as t~ 3 / 2 . It should cross the 8.46 GHz wavelength around day 140. This is much later 
than the time to for transition to slow cooling, which occurs around to = ^.le^^A* days, 
for typical parameters. Therefore, we can use eq. (2) to find the expected peak flux at 8.46 
GHz. The result is 

^ m ,8.46GHz = O.SQe^iiA, mJy, (4) 

which for standard parameters is comparable to the R-band flux at the break t m ^. The 
8.46 GHz flux is observed at 598±33 yuJy on day 5.7. If this flux is comparable to the 
peak flux at the time t m (which is true if the observed frequency is close to or beyond the 
self- absorption frequency; see CL00), then one can use eq. (4) to provide a rough estimate 
for A* in terms of e#: A* m l.Oe^ 2 / 3 . 

The above analytic estimates assumed p = 2.5, and are rather crude. They do indicate 
that the optical and radio data may be fitted with a wind model with parameters not far 
from the standard values. We now demonstrate that this is indeed the case using a numerical 
model. The model treats synchrotron emission from a spherical (trans-)relativistic blast 
wave propagating in an r~ 2 density medium, with the light travel time effects, synchrotron 
self absorption and cooling included. It was previously applied to GRB 980508 (CL00), 
among others. Like GRB 980508, the decay of the R-band flux on the time scale of days 
and longer is relatively slow, with F u oc t^ 1 approximately; the wiggles on the light curve 
of this source makes a precise determination of the decay slope difficult. The slope implies 
that p is close to 2, although the exact value is somewhat uncertain. We pick p — 2.1, which 
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corresponds to a decay slope of a = —(3p — 2)/4 = 1.075 in the optical. 

After some experimentation, we find a solution that fits the R-band and radio data 
reasonably well with the following set of parameters: e e = 0.1, e B = 0.1, E 52 = 10, and 
A* = 0.6, which corresponds to a wind mass loss rate of 6 x 10~ 6 M yr^ 1 (assuming a 
nominal wind speed of 10 3 km s _1 ). Now all the parameters are approximately determined 
(within a factor ~ 2) because of the inclusion of self-absorption. The fits are shown in 
panels (a) and (b) of Fig. 1. We did not attempt to fit the bumps on the R-band light curve; 
they have been interpreted as arising from either late energy injections or inhomogeneities 
in the ambient density (Lazzati et al. 2002; Nakar, Piran & Granot 2003; Heyl & Perna 
2003). The bumps introduce some uncertainty to the model parameters we obtained. Radio 
emission was detected at 4.86, 8.46, 15, 22.5 and 86 GHz at various times (Frail & Berger 
2002; Berger, Frail & Kulkarni 2002; Pooley 2002b,c; Bermer & Castro-Tirado 2002). One 
upper limit each exists at 15, 232 and 347 GHz (Pooley 2002a; Bermer & Castro-Tirado 
2002; Wouterloot et al. 2002). The flux measurement of 2.5±0.3 mJy at 86 GHz flux at an 
average time of 1.5 days is particularly interesting. This flux is three times higher than the 
R-band flux at the break around day 0.1. It presents a problem to the identification of the 
break as t m in a constant density medium when the cooling frequency has already passed 
the R-band from below (Kobayashi & Zhang 2002). In such a case, the maximum fluxes 
at lower frequencies should be the same as that of the R-band at the break, namely about 
0.83 mJy, which is well below the 86 GHz measurement. This discrepancy was also noted 
by Pandey et al. (2002). 

The relatively high 86 GHz flux is not a problem for our model, where the transition 
to slow cooling occurs around day 2, much later than day 0.1. In a wind model, the flux 
can be much higher in the radio (broadly defined to include millimeter and sub-millimeter 
wavelengths), particularly at early times when the cooling frequency is expected to be in 
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the spectral region. This behavior shows up clearly in panel (b) of Fig. 1, where the peak 
fluxes at the three highest frequencies are all above 3 mJy. Such high fluxes are naturally 
expected in a wind model but not in an ISM model, as emphasized by Panaitescu & Kumar 
(2000). However, it is difficult to make a strong case for wind interaction based on a single 
data point at 86 GHz. A stronger case can be made if the 8.46 GHz flux starts to decline 
around 100 days, when the typical frequency v m is expected to pass through the frequency 
from above. This expectation needs to be modified in the case of an early jet break. 

X-ray afterglows are observed with Chandra at two epochs. The first epoch started 
about 0.87 days after the burst and lasted for 88.1 ksec (Sako & Harrison 2002a; Fox et 
al. 2003b). Within this epoch, the X-ray afterglow has a power-law spectrum, with index 
Px = —1-1 ± 0.1, and decays roughly as a power-law, with index ax = —1.0 ± 0.2. Both are 
consistent with our model, where the cooling frequency u c around day 1 is well below the 
X-ray band, and fix and «i are predicted to be -1.05 and -1.075, respectively, for p = 2.1. 
The fact that the temporal decay slope in the R-band, «o, is close to -1 indicates that the 
cooling frequency is below the R-band as well around this time, which is in agreement with 
the optical spectral index of Po = —1.07 ± 0.06 determined by Pandey et al. (2002) and 
Po — —0.96 ± 0.03 by Matheson et al. (2003) in the absence of a substantial host galaxy 
extinction (see also Bersier et al. 2003; Holland et al. 2002). Our best fit model yields an 
X-ray flux of 2.6 x 10~ 13 erg cm" 2 s" 1 between 2 and 10 keV at day 1.4, close to the middle 
of the first observing epoch. It is lower than, but within a factor of two of, the mean flux of 
the entire epoch 4.3 x 10~ 13 erg cm" 2 s" 1 (Sako & Harrison 2002a). We therefore conclude 
that the wind model is consistent with the first epoch of X-ray observations. 

The second epoch of X-ray observations started 52.23 days after the burst and yielded 
a 2-10 keV flux of 7.2±2.5x 10~ 16 erg cm" 2 s" 1 (Sako & Harrison 2002b). The flux implies a 
decay slope between the two epochs of approximately ax = —1-7, which is steeper than that 
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predicted in our spherical model. The steepening may be due to a jet break between the 
two epochs of observation, which was suggested to have occurred around day 7 by Pandey 
et al. (2002) and Holland et al. (2002) based on their interpretation of the (wiggly) R-band 
light curve. The jet break, if exists, should show up in a well-sampled radio light curve as 
well. Alternatively, the steepening could be due to a steepening of the energy distribution 
of nonthermal electrons well above the minimum energy of the electrons accelerated at the 
shock front (e.g., Li & Chevalier 2001; Panaitescu & Kumar 2002). 

The question of a jet break is related to the energy in the source. The energy we find 
in a spherical model for the afterglow, E 52 = 10, is comparable to the isotropic burst energy 
in 7-rays, 5 x 10 52 ergs (Bloom, Frail, & Kulkarni 2003). Pandey et al. (2002) find, in a fit 
to the optical data through day 21, that there is a break in the light curve at t b = 7.6 ± 0.3 
days. Bloom et al. (2003) interpret this as a jet break; the correction for collimated flow 
reduces the 7-ray energy by a factor of ~ 40. However, Fig. 1 shows that the variability 
in the light curve and the late light curve points make a clear designation of the jet break 
difficult. In a wind model, the jet break evolves slowly (Kumar & Panaitescu 2000), which 
makes any determination of a jet break from afterglow data more uncertain. 

The mass loss rate we deduced is typical for a Wolf-Rayet type wind. There are other 
indications that GRB 021004 may be interacting with a Wolf-Rayet type wind. Wolf-Rayet 
winds are thought to be clumpy (e.g., Hamann & Koesterke 1998), and the dumpiness 
may provide an explanation for the prominent bumps on the R-band light curve of GRB 
021004 (Lazzati et al. 2002; Nakar et al. 2002; Heyl & Perna 2003). In addition, there 
are multiple absorption components in the spectrum of the afterglow, separated by speeds 
up to 3000 km s -1 . These components could come from substructures in a Wolf-Rayet 
wind (Mirabal et al. 2002; Schaefer et al. 2002). Interestingly, GRB 990510, which is best 
modeled by interaction with a constant density medium (CL00; Panaitescu & Kumar 2002), 
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shows a smoothly evolving optical afterglow (Stanek et al. 1999). 

GRB 021004 differs from the other two GRBs with detected early afterglows (GRBs 
990123 and 021211) in several ways: it has a slow decay in the R-band light curve followed 
by a steepening rather than a steep decline followed by a flattening (which occurs at a 
much earlier time than the break in GRB 021004), a higher optical flux at late times after 
the break, and a bright, long-lived radio afterglow. The early emission from GRB 990123 
was convincingly interpreted as coming from the reverse shock of a blast wave expanding 
into an ISM (Sari & Piran 1999; Meszaros & Rees 1999), and the late-time afterglow data 
are consistent with an ISM model (e.g., Panaitescu & Kumar 2002). The afterglow of GRB 
021211 resembles that of GRB 990123, and was interpreted similarly (Fox et al. 2003a; Li 
et al. 2003; Wei 2003). The rate of initial decline in the R-band flux of approximately t~ 2 
(GRB 990123) or shallower (GRB 021211) is difficult to reproduce in the reverse shock of 
a wind model. The available data on the three early afterglows therefore appear to point 
to two types of GRB ambient environments: ISM and a stellar wind. These two types of 
environments are also inferred in the detailed modeling of a larger number of afterglows at 
later times (CL00; Panaitescu & Kumar 2002; Price et al. 2002). 

The evidence for wind and constant density environments in the early afterglow 
evolution supports the hypothesis that these represent different progenitors (Chevalier 
& Li 1999; CL00) because a small radial distance from the star is explored (^ 10 17 cm) 
where a freely expanding wind is expected. A prediction of the CL99 scenario is that wind 
interaction should be correlated with supernova light; this hypothesis was supported by 
GRB 011121 (Price et al. 2002). However, GRB 021211, which has an early afterglow 
indicating ISM interaction, shows some evidence for a supernova-like bump in the light 
curve (Fruchter et al. 2002; Testa et al. 2003). GRB 020405 is another case of apparent 
ISM interaction and a supernova-like bump (Berger et al. 2003). Spectra confirming the 
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supernova nature of light curve bumps, as well as detailed afterglow observations, are 
needed to clarify the situation. Long term monitoring of radio afterglows will be crucial 
in testing the wind model of early afterglows (by examining the evolution of v m ), and in 
determining the wind density. Another prediction of the wind model is that the early 
optical emission before the break (t < t m ) should have the spectrum F u oc z/ -1 / 2 , which is 
a flatter spectrum than is typically observed in optical afterglows. In addition, the break 
is chromatic, occurring at a later time for a longer wavelength, which can be tested with 
densely sampled IR observations. 

Support for this work was provided in part by NASA. 
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time since explosion (days) 



Fig. 1. — Wind interaction model for the afterglow of GRB 021004. The optical data are 
taken from the papers Fox et al. (2003b), Bersier et al. (2003), and Holland et al. (2003) and 
the GCN notices Matsumoto et al. (2002), Weidinger et al. (2002), Mirabal et al. (2002) 
and Fatkhullin et al. (2002), with the modest amount of Galactic extinction corrected. The 
radio data are taken from Frail & Berger (2002; 22.5 GHz), Berger, Frail & Kulkarni (2002; 
4.86 and 8.46 GHz), Pooley (2002b,c; 15 GHz), and Bermer & Castro-Tirado (2002; 86 
GHz). The upper limits at 15, 232 and 347 GHz are given in Pooley (2002a), Bermer & 
Castro-Tirado (2002) and Wouterloot et al. (2002). The lines are the light curves from the 
wind model described in the text. 



